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The  contract  has  supported  an  extensive  programme  of  atmospheric 
and  ionospheric  modelling,  aimed  at  ezLaminlng-^tmospb.q,r_ij:  Structure  and 
Variability,  using  the  CRAY  XMP-48  "at  the  Rutherford  Applet^rTLab'ora'Ctiryj-nh'd '  the 
CRAY  2S  machines^at  the  University  of  London  Computing  Centre,, part  of  the  ' 
programme  o^atffiospherlc  physics  research  activities  of  the  Atmospheric  .Physics 
Laboratory.  A  technique  of  forcing  tidal  and  other  natural  atmospheric 
tluctuatlons  at  the  lower  boundary  of  the  present  version  of  the  3- dimensional 
time-dependent  global  thermosphere -ionosphere  model  has  been  developed.  Also,  a 
series  of  numerical. studies  have  been  carried  out  into  the  modelled  response  of 
the  thermosphere  during  a  series  of  geomagnetic  storms  in  1981  and  1982,  observed 
by  the  NASA  Dynamics  Explorer-2  spacecraft.  Further  studies  have  concentrated  on 
using  these  models  for  studies  of  thermosphere  -  ionosphere  coupling,  atmospheric 
density,  lower,  atmospheric  tides,  participation  in  the  analysis  and  interpretation 
Of  empirical  data  obtained  during  the  'Lower  Thermosphere  Coupling  Study',  the 
mid-  and  low-latitude  ionospheric  effects  of  geomagnetic  storms,  and  for  studies 
of  the  solar  activity  and  geomagnetic  activity  causes  of  variations  of  lower 
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SUMMARY 

An  extensive  programme  of  atmospheric  and  ionospheric  modelling,  aimed  at  examining 
"Atmospheric  Structure  and  Variability",  using  the  CRAY  XMP-48  at  the  Rutherford 
Appleton  Laboratory  and  the  CRAY  2S  machines  at  the  University  of  London  Computing 
Centre  has  been  carried  out  as  part  of  the  atmospheric  physics  research  programme 
of  the  Atmospheric  Physics  Laboratory.  A  summary  of  our  major  projects, 
highlighting  the  recent  research  achievements  and  publications,  is  given  below. 
There  have  been  two  major  thrusts  of  the  .ork  in  progress  under  this  contract. 

Firstly,  the  development  of  a  technique  of  forcing  tidal  and  other  natural 
atmospheric  fluctuations  at  the  lower  boundary  of  the  present  version  of  the  3- 
dimenslo,ial,  time-dependent  global  thermosphere -ionosphere  model.  This  Improves  the 
simulation  of  the  behaviour  of  structure  and  variability  within  the  upper  mesosphere 
and  lower  thermosphere.  With  this  tidal  model  It  is 'possible  to  Intercompare  the 
model  results  with  data  from  ground-based  instruments,  such  as  the  incoherent 
scatter  and  K.S.  radars  and  Fabry-Perot  Interferometers. 

Secondly,  a  series  of  numerical  studies  have  been  carried  out  into  the  modelled 
response  of  the  thermosphere  during  a  series  of  geomagnetic  storms  in  1981  and  1982. 
These  storms  were  well-observed  by  the  NASA  Dynamics  Explorer-2  spacecraft,  and  it 
his  been  possible  to  carefully  compare  the  series  of  model  simulations,  using 
various  parameterised  drivers  representing  geomagnetic  variability.  The  response  of 
the  upper  thermosphere  to  geomagnetic  activity  is  one  of  the  major  causes  of  density 
temperature  and  wind  variability.  The  most  appropriate  index  of  activity  is  quite 
central  to  improving  the  modelling,  using  empirical  or  theoretical  models.  We  have 
shown  that  the  difference  between  the  time-dependence  of  available  indices  is  quite 
critical  for  accurate  modelling,  so  that  there  is  as  large  a  difference  between 
model  results  using  different  available  Indices,  as  between  any  of  the  storm 
simulations  and  data  from  Dynamics  Explorer.  Investigation  of  the  causal  physics 
via  the  combination  of  observations  (l.e.  DE,  and  UARS  and  other  future  missions) 
and  simulations  using  the  present  coupled  ionosphere-thermosphere  models,  will 
improve  our  understanding  of  how  to  best  combine  the  available  information  on 
thermosphere  /  ionosphere  drivers.  We  expect  this  to  lead  to  Improved  indexing  of 
both  the  empirical  and  theoretical  models  to  geomagnetic  activity. 

The  large  "time  allocations  made  available  on  the  UK  CRAY  machines  have  been 
effectively  used  to  Implement  our  theoretical  and  numerical  modelling  research 
programme  during  the  past  2  years .  This  work  has  been  primarily  concerned  with  the 
extended  development  of  the  UCL  /  Sheffield  coupled  model  of  the  neutral 
thermosphere  and  ionosphere.  This  has  Involved  the  integration  of  the  University 


College  London  neutral  atmosphere  code,  designed  to  model  the  energy  budget, 
dynamics,  and  composition  of  the  earth’s  upper  atmosphere,  with  the  Sheffield 
University  high- latitude  Ionospheric  model.  The  latter  model  Includes  the  effects 
of  magnetospherlc  plasma  convection,  energetic  particle  precipitation  and 
photolonlsatlon.  The  coupled  model  Is  able  to  self -consistently  simulate  the 
variety  of  complex  Interactions  between  the  neutral  and  plasma  environments .  This 
enables  us  to  examine  the  importance  and  relevance  of  the  different  processes,  and 
to  evaluate  their  individual  roles  In  producing  observed  phenomena.  The  Ionospheric 
code  has  recently  been  extended  to  cover  mid  and  low  latitude  regions.  We  have  run 
the  combined  model  on  the  CRAY  machines  and  have  made  a  number  of  significant 
advances  In  the  subject  during  the  past  2-3  years. 

RECENT  PROGRESS 

Significant  progress  his  been  made  In  the  two  major  aspects  of  the  modelling  work: 

Firstly,  there  are  the  Interactions  with  tidal,  planetary  and  gravity  waves 
emanating  from  the  lower  and  middle  atmosphere,  which  primarily  affect  the  lower 
boundary  conditions  and  the  structure  and  dynamics  of  the  lower  thermosphere.  With 
the  current  model,  we  can  now  intercompare  rest  Its  from  the  incoherent  scatter 
radars  and  new  generation  ground-based  FPIs,  this  has  lead  to  a  number  of  papers  now 
In  press  and  in  preparation.  The  empirical  studies  show,  very  convincingly,  the  day 
to  day  variability  of  the  tidal  and  gravity  wave  forcing  of  the  lower  thermosphere 
via  troposphere  /  stratosphere  sources,  providing  a  new  challenge  for  the  extended 
development  of  the  numeric  il  codes. 

Secondly,  we  have  continued  studying  the  Interactions  between  the  magnetosphere  and 
the  polar  Ionosphere,  through  the  medium  of  the  polar  electric  field,  and  the 
magnetospherlc  field-aligned  current  system,  with  charge  carriers  Including  the 
energetic  particles  which  excite  the  auroral  Ionosphere  and  auroral  emissions.  The 
first  experiments  with  the  Rice  University  Group  have  now  been  performed. 
Interchanging  model  results,  so  chat  we  can  examine  the  nature  of  Interactions, 
ahead  of  fully -coup ling  the  two  numerical  codes.  Even  the  present  computing 
machines  have  a  critical  limitation  In  the  core  and  speed  requirements  for  such  a 
task.  However,  new  machines,  expected  to  be  available  for  'pump -priming’  by  mid 
1990,  will  be  able  to  run  such  complex  tasks  as  a  fully -Interactive  magnetosphere- 
lonosphere-thermosphere  model.  We  anticipate  that  this  activity  will  lead  to 
Improved  understanding  of  the  relationship  between  the  actual  magnetospherlc  sources 
of  energy  and  momentum,  their  coupling  to,  and  sinks  within  the  thermosphere,  and 
the  various  feed-back  processes.  We  also  hope  that  this  work  will  lead  to  a  deeper 
understanding  of  the  physical  processes  behind  the  available  Indices  of  global 
geomagnetic  activity  and  solar  wind  parameters,  which  will  aid  future  predictive 
modelling  of  the  Ionosphere -thermosphere  system  and  its  response  to  solar  and 
geomagnetic  variability. 

I .  Thermosphere/Ionosohere  coupling. 

Simulations,  by  the  coupled  thermosphere/ionosphere  model,  for  the  December  and  June 
solstices  have  shown  perhaps  the  most  striking  Interactions  between  the  neutral  and 
Ionised  components .  The  model  has  shown  that  In  the  sujamer  polar  regions ,  the  upper 
thermosphere  neutral  gas  composition  has  a  very  strong  control  on  the  plasma  density 
/!/.  This  strong  control  was  not  apparent  when  using  recent  empirical  models  of 
neutral  composition  such  as  MSIS  1986.  Comparison  of  the  seasonal  variation  of 
plasma  density,  as  observed  by  the  European  Incoherent  Scatter  Radar  facility 
(EISCAT)  In  northern  Scandinavia,  now  shows  good  agreement  with  the  simulations  /2/. 
The  model  now  enables  us  to  understand  relative  Importance  of  the  different  major 
processes  which  are  Involved. 


I 

The  simulations  of  the  electrodynamics  of  the  lower  thermosphere  at  high  latitudes 
are  now  realistic,  enabling  the  model  to  be  used  to  test  the  various  empirical 
descriptions  of  the  magnetospherlc  energy,  momentum  and  precipitation  inputs  /3/, 
which  are  some  of  the  main  drivers  of  the  upper  atmosphere.  In  addition,  the 
simulated  neutral  circulation  In  the  lower  thermosphere  at  high  latitudes  In  these 
self-consistent  calculations  Is  much  closer  to  observations  /4,5/. 

During  equinox  and  at  winter  solstice,  transport  of  plasma  by  the  convective 
electric  field  has  a  major  Influence  on  the  structure  of  the  high  latitude 
Ionosphere.  The  nvimerlcal  model  has  been  used  to  Isolate  the  Influence  and 
respective  contributions  of  the  polar  convective  electric  field  and  auroral 
precipitation  In  determining  the  structure  of  the  ionospheric  F-region,  and  In  the 
creation  of  the  trough  fea^ureu  oi>  the  equacor.a.d  edge  of  the  auroral  oval  /6/. 


II .  Atmospheric  Density. 

The  coupled  numerical  model  provides  a  more  accurate  means  of  computing  the  response 
of  the  density  of  the  upper  atmosphere  to  changes  In  solar  radiation,  the 
magnetospherlc  Input,  and  from  tides  propagating  from  the  lower  thermosphere.  An 
accurate  knowledge  of  the  mass  density  Is  necessary  for  accurate  prediction  of  the 
location  and  lifetimes  of  earth  orbiting  satellites.  Simulations  have  been 
performed  comparing  the  various  available  Indices  of  geomagnetic  and  auroral 
activity.  These  Indices  describe  the  magnetospherlc  input  in  three  distinct  ways 
/7/.  The  results  Indicate  that  prediction  of  small  scale  and  short  period 
variations  requires  considerably  Improved  empirical  data  to  describe  the  geomagnetic 
Inputs  and  their  variations. 

ill.  Atmospheric  Tides . 

A  major  driver  of  the  mid  and  low  latitude  thermosphere  is  the  tidal  energy 
propagating  upwards  from  sources  In  the  lower  atmosphere.  The  thermosphere  is  slow 
to  respond  to  this  forcing  (5  to  10  days),  requiring  long,  time-consuming, 
computation.  Numerical  simulations  of  the  seasonal  response  to  propagating  tides 
Indicates  that  the  magnitudes  and  phases  of  propagating  tides  at  the  base  of  the 
thermosphere  are  quite  variable.  This  variable  forcing  is  required  to  match 
empirical  data  obtained  from  different  locations  and  observing  periods  /8,9/.  We 
have  recently  demonstrated  that,  with  our  new  Fabry-Perot  interferometer  Installed 
at  the  Utah  State  University  Hardware  Ranch  Facility  at  Bear  Lake,  Utah,  we  can 
obtain  very  high  quality  data  on  tidal  wind  components  around  86  km  altitude,  using 
the  Helnel  band  emissions  at  843  nm  /lO/.  These  continuing  tidal  wind  studies, 
combi.  Ing  new  empirical  and  model  studies,  are  described  further  In  Annex  1. 

IV.  Lower  Thermosphere  Coupling  Study 

Comparisons  have  been  made  of  the  model  prediction  of  the  lower  thermosphere  neutral 
wind  and  temperature  and  Ion  temperature  with  Incoherent  Scatter  observations  for 
the  September  1987  Lover  Thermosphere  Coupling  Study  /9/. 

The  results  Indicate  that  symmetric  (2,2)  and  (2,4)  global  Hough  modes  can  reproduce 
the  vertical  structure  of  the  observations  at  high  and  mid-latitude  provided 
amplitude  significantly  larger  than  linear  model  predictions  are  used.  The 
simulations  also  show  that  at  high  latitudes  the  seml-dlumal  response  In  Ion 
temperature  Is  driven  primarily  by  geomagnetic  processes,  rather  than  tides,  and 
that  a  large  difference  In  the  Ion  and  neutral  temperature  penetrates  down  to  110  km 
altitude  during  active  conditions. 
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V.  Geomagnetic  Storm  Effects 


Observations  from  mid-latitude  lonosondes  indicate  that  the  F2  layer  electron 
density  Is  depleted  during  and  following  a  geomagnetic  storm.  The  'negative'  storm 
phenomenon  is  a  complex  and  Important  consequence  of  geomagnetic  storms  and  affects 
global  radio  communications.  The  coupled  model  has  been  used  to  assess  the 
consequences  of  a  possible  storm-time  heat  source  resulting  from  energetic  ion 
precipitation  from  the  ring  current.  It  was  shown  that  such  an  additional  source 
could  appreciably  change  the  neutral  composition,  enhancing  the  F2  region  depletion 
/ll/. 

Recent  simulations  with  the  coupled  thermosphere  Ionosphere  model  have  been  able  to 
reproduce  the  seasonal  and  local  time  response  of  thermospheric  composition  to 
storm-time  magnetospheric  Input  /12/.  The  storm-time  model  simulations  reproduce 
the  penetration  to  low  latitudes  of  a  region  of  increased  molecular  nitrogen  density 
in  the  early  morning  sector  of  the  summer  hemisphere.  These  changes  qualitatively 
reproduce  the  observed  depletions  in  the  F2-region  Ionosphere  expected  as  the  result 
of  an  Increase  In  mean  molecular  mass. 

VI .  Atomic  oxygen  and  nitric  oxide. 

Using  an  advanced  version  of  a  2-Dlmensional  global  model,  which  includes  a  more 
comprehensive  photo-chemical  code  than  the  3-D  model,  it  has  been  possible  to  study 
the  global  distribution  of  atomic  oxygen  and  nitric  oxide  as  a  function  of  season, 
latitude  and  geomagnetic  activity  /13/.  The  data  recently  published  on  nitric  oxide 
distributions  and  variability  from  SME  satellite  observations  by  Barth  (1988)  and 
Siskind  et  al  (1989)  has  provided  an  excellent  test  of  the  model,  and  an  opportunity 
to  use  the  model  to  test  the  key  branching  ratios  involved  in  the  photochemical 
production  processes,  and  the  wide  range  of  atmospheric  processes  involved  In  the 
generation,  distribution  and  destruction  of  nitric  oxide,  and  its  implications  for 
the  distribution  of  other  key  minor  constituents,  including  mesospheric  ozone. 

VII.  Review  papers. 

Results  from  the  studies  carried  out  under  the  support  from  this  contract  have  been 
used  in  a  number  of  other  invited  review  papers  presented  in  recent  International 
conferences  (AGARD,  Chapman  'Auroral  Physics'  Conference  in  Cambridge,  July  1988, 
COSPAR  (July  1988,  NATO  (Sept  1988,  Moscow  (Jan  1989),  ESA  /ESTEC  (May  1989)  and 
lAGA,  Exeter,  July  1989,  /12-22/. 

VIII.  Other  Work. 

The  UCL  /  Sheffield  Model  has  been  used  In  conjunction  with  the  interpretation  of 
data  obtained  during  a  number  of  coordinated  ground-based  programmes  within  WITS, 
such  at  GITCAD,  LTCS ,  ETC.  Approximately  10  papers  have  been  recently  published  or 
are  in  press  describing  the  results  of  these  joint  programmes  and  studies. 
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ANNEX  1. 

SUMMARY  OF  RELEVANT  WORK  ON  MIDDLE  AND  UPPER  ATMOSPHERE  TIDES  USING  THE  UCL  MODEL 


Coupling  of  tides  generated  within  the  troposphere,  stratosphere  and  mesosphere  and 
interacting  with  the  in-situ  solar  tides  and  the  results  of  geomagnetic  forcing  of 
Che  thermosphere  has  been  simulated  within  the  UCL  3-diraenslonal  time -dependent 
model  of  the  thermosphere  via  the  introduction  of  tidal  forcing  near  the  lower 
boundary  of  the  existing  coupled  ionosphere  -  thermosphere  model  (80  km)  . 

The  present  UCL  coupled  model  thus  has  been  used  to  simulate  Che  amplitudes,  phases 
and  height  structure  of  the  combined  tides  in  the  thermosphere.  With  appropriate 
tuning  of  the  tidal  sources,  results  obtained  are  in  reasonable  agreement  with  tidal 
wind  observations  at  low,  mid  and  high  latitudes. 


The  influence  of  propagating  tides  on  the  thermosphere  has  been  studied  for  varying 
background  conditions,  as  a  function  of  latitude,  season,  solar  activity  and 
geomagnetic  activity.  The  significant  interactions  between  propagating  and  in- situ 
thermospheric  tides  and  between  different  propagating  tidal  modes  have  also  been 
investigated.  Some  of  these  interactions  are  because  of  the  inherent  non-linear 
coupling  of  various  tidal  modes  within  the  mesosphere  and  lower  thermosphere.  under 
other  conditions,  it  is  the  presence  of  the  ionosphere  which  leads  to  non-linear 
behaviour  and  tidal  mode  coupling. 

We  have  recognised  from  the  difficulty  of  fitting  widely  different  sets  of  tidal 
wind  data  with  a  single  description  of  the  source  (represented  by  a  geopotential 
variation  introduced  near  90  km  altitude)  that  large  variations  in  the  propagation 
conditions  and  characteristics  within  the  stratosphere  and  mesosphere  probably 
occur.  Such  changes  can  affect  considerably  the  amplitudes  and  phases  of  the 
propagating  tides  well  before  they  reach  the  lower  thermosphere.  The  effects  ot 
such  changes  on  the  lower  thermospheric  winds  have  been  simulated  by  varying  the 
input  amplitude  and  phase  ot  tidal  forcing  near  the  existing  lower  boundary  of  the 
thermosphere  model. 

All  of  these  simulations  of  propagating  tides  have  used,  for  the  first  time,  a 
fully-coupled  ionosphere  -  thermosphere  model.  In  this  model,  all  of  the  raoraentura, 
chemical  and  dynamical  interactions  between  the  ionised  and  neutral  components  are 
included  have  been  modelled  in  a  self-consistent  manner,  something  which  has  not 
been  possible  in  previously  modelling  attempts  (Forbes  and  NCAR) .  Naturally,  the 
momentum  interchange  effects  are  strongest  in  the  geomagnetic  polar  regions, 
however,  important  effects  from  the  point  of  view  of  tidal  coupling  occur  at  all 
latitudes , 

The  relative  importance  of  geomagnetic  and  tidal  forcing  mechanisms  at  high 
latitudes  has  already  been  investigated  using  this  coupled  model.  The  results  of 
model  simulations  at  high  latitude  have  been  compared  with  observations  which  have 
only  recently  become  available  from  the  EISCAT  incoherent  scatter  radar  facility. 

The  (1,1)  diurnal  propagating  tidal  mode,  which  is  important  in  the  lower 
thermosphere,  has  also  been  simulated  for  the  first  time  in  any  three-dimensional 
time-dependent  thermospheric  model.  This  work  has,  of  course,  used  the  coupled 
ionosphere  -  thermosphere  model,  which  is  also  a  first.  The  existence  of  interactions 
between  the  diurnal  and  semidiurnal  tides  have  also  been  demonstrated  within  the  UCL 
mode  1 

The  interactions  of  propagating  tides  with  the  zonally-avoraged  flow  have  also  been 
studied.  The  dissipation  of  propagating  tides  has  been  shown  to  produce  strong 
zonal  accelerations.  Such  zonal  accelerations  have  previously  been  found  in  models 
of  gravity  wave  dissipation.  The  relative  importance  of  different  dissipation 
mechanisms  for  propagating  tides  as  a  function  of  height  have  also  been  investigated 
within  the  UCL  model. 

Recently,  we  have  obtained  some  fascinating  data  on  tides  within  the  mesosphere  from 
the  new  Fabry-Perot  interferometer  observing  hydroxyl  Meinel  band  emissions  at  Utah 
State  University,  These  data,  which  show  very  large  (50  ms-1)  semidiurnal  tides, 
some  excellent  examples  of  gravity  waves,  and  occasional  very  large  wind 
disturbances,  provide  many  challenges. 


OUTLINE  PROPOSAL  FOR  FURTHER  DEVELOPMENT: 


In  response  to  these  challenges,  we  are  currently  seeking  the  resources  for 
additional  development  of  the  coupled  ionosphere  -  thermosphere  model: 

i) .  To  decrease  the  altitude  of  the  lower  boundary  of  the  coupled  model  to  40  kin; 

ii)  To  introduce  self-consistent  geopotential  variations  just  above  the  40  kra  lower 
boundary  to  simulate  self -consistent  tidal  forcing  of  the  entire  mesosphere  and 
thermosphere  (above  60  -  65  km); 

iii)  To  develop  a  comparable  scheme  for  self-consistent  introduction  of  gravity 
waves  and  other  disturbances  comparable  Co  chose  we  are  currently  observing  with  the 
UCL  /  USU  FPI  at  Hardware  Ranch; 

iv)  To  introduce  self-consistent  geopotential  forcing  near  the  new  lower  boundary 
to  simulate  the  effects  of  Planetary  wave  structures  of  the  mesosphere. 

It  is  clear  that  the  tidal  and  gravity  waves  which  are  observed  cannot  be  adequately 
simulated  by  assuming  that  the  propagation  characteristics  within  the  upper 
stratosphere  and  mesosphere  are  always  unchanged.  Thus  there  is  the  challenge-  to 
handle  that  variability  in  a  fully  self-consistent  fashion.  The  fabric  and 
structure  of  the  UCL  coupled  model  is  designed  for  such  adaptation,  and  we  art- 
seeking  the  additional  manpower  resources  to  aid  us  accomplish  Chat  task, 
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